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Abstract. We study the transition in conductance properties of chaotic mesoscopic 
cavities as time-reversal symmetry is broken. We consider the Brownian motion model 
for transmission eigenvalues for both types of transitions, viz., orthogonal- unitary and 
symplectic-unitary crossovers depending on the presence or absence of spin-rotation 
symmetry of the electron. In both cases the crossover is governed by a Brownian 
motion parameter r, which measures the extent of time-reversal symmetry breaking. 
It is shown that the results obtained correspond to the Jacobi crossover ensembles 
of random matrices. We derive the level density and the correlation functions of 
higher orders for the transmission eigenvalues. We also obtain the exact expressions 
for the average conductance, average shot-noise power and variance of conductance, as 
functions of r, for arbitrary number of modes (channels) in the two leads connected to 
the cavity. Moreover, we give the asymptotic result for the variance of shot-noise power 
for both the crossovers, the exact results being too long. In the r — ► and r — > oo 
limits the known results for the orthogonal (or symplectic) and unitary ensembles are 
reproduced. In the weak time-reversal symmetry breaking regime our results are shown 
to be in agreement with the semiclassical predictions. 
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1. Introduction 

The transport of electrons through ballistic mesoscopic cavities has emerged as a subject 
of intensive research in recent years [U[2l[3l[ll[5l[6l[Tl[8l[9l HOI HH H21 [131 HS1 HH1 [TTl 
UHl Uni EQl EH [221 [231 EH EHJ [26l EZl EH [29l EQl EU [32]- The technological advancement 
has made it possible to fabricate mesoscopic cavities of arbitrary geometries and degree 
of disorder (TJ |2j [3j HJ E]. The main theoretical tools used to study the quantum 
conductance phenomenon in these systems include random matrix theory (RMT) [3], 
diagrammatic methods [7J , semiclassical methods (8] and supersymmetry techniques [9] . 
RMT provides a convenient way to study the transport properties when the underlying 
classical dynamics of these cavities is chaotic [U EJ [3j HJ [TOj E] , and it has been quite 
successful. 

The scattering matrix S, which relates the electron fluxes in the two leads connected 
to the mesoscopic cavity, is central to the problem of quantum transport. Under RMT 
the scattering matrix is assumed to belong to one of Dyson's standard circular ensembles. 
A microscopic justification for using circular ensembles as models for scattering matrices 
has been provided by calculations involving the Gaussian ensembles [T2l H3] . The 
appropriate circular random matrix ensemble for the problem is decided by the time- 
reversal and rotational symmetries of the system and is determined by its invariance 
under orthogonal, symplectic or unitary transformations [3] [33]. The invariant ensembles 
are abbreviated as OE, SE and UE respectively and characterize the three classes of 
chaotic cavities. 

Various physical quantities of interest such as average conductance, average shot- 
noise power and variance of conductance derive from the statistics of transmission 
eigenvalues, which are obtained by considering the polar decomposition of the scattering 
matrix. It turns out that the statistics of transmission eigenvalues Tj for the invariant 
cases is related to the statistics of eigenvalues of a special case of Jacobi ensembles under 
a simple linear transformation of Tj [3j |30j |31j [32]. The relation between transmission 
eigenvalues obtained from the scattering matrix and the eigenvalues of Jacobi ensembles 
of random matrices has been explicitly shown in [32J for the invariant cases. 

The problem of statistics intermediate between above invariant classes is also 
important. For instance, application of a weak magnetic field to the system breaks the 
time-reversal symmetry partially. In this case, depending on the presence or absence of 
spin-rotation symmetry of the electron, statistics intermediate between OE and UE or 
SE and UE is observed. The problem of intermediate statistics in OE-UE and SE-UE 
transitions has been solved exactly for the energy level correlations as modeled by the 
Gaussian [Ml [35] and circular [36] ensembles and recently for the entire Jacobi family 
of random matrix ensembles [37] . 

Our purpose in this paper is to study the intermediate statistics for the conductance 
in mesoscopic systems. This has been a subject of interest both theoretically and 
experimentally [201 EH E21 E3 EH EH1 EH [27l EHl ES]- The problem of intermediate 
statistics is studied under the framework of Dyson's Brownian motion model of 
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random matrices [351 EH] • The corresponding Brownian motion model for transmission 
eigenvalues has been obtained by Frahm and Pichard [20J using the polar decomposition 
of the scattering matrices. This problem of finding the transmission eigenvalue statistics 
for the crossovers is technically challenging, because there is no direct relation between 
the transmission eigenvalues and the eigenphases of scattering matrix [36]. Frahm and 
Pichard have considered the model for equal number of incoming and outgoing channels 
and solved the OE-UE crossover problem [21]. We extend their formalism to unequal 
number of incoming and outgoing channels for both OE-UE and SE-UE crossovers. We 
show that the results for the transmission correlation functions correspond to Jacobi 
crossover ensembles [37]. We use these results to calculate the above-mentioned physical 
observables as functions of the Brownian motion parameter r. In the r — > and r — > oo 
limits, results corresponding to OE (or SE) and UE, are reproduced. Also in the weak 
time-reversal symmetry breaking regime our results are consistent with the semiclassical 
predictions [25, 26]. 

The paper is organized as follows. In section [2] we make use of the Brownian 
motion model to obtain the jpd of transmission eigenvalues in the crossover regimes. 
In section [3] we introduce the two-point kernels which are used to calculate the n-level 
correlation functions. Also, we express the averages and variances of conductance and 
shot-noise power in terms of moments involving these kernels. In section H] we consider 
the OE-UE crossover and give the explicit expression for the density of transmission 
eigenvalues. This generalizes the result available for the UE [301 EI] to the crossover. 
Using the Landauer-Biittiker formalism [HI [151 E3 EE] we find the exact results for 
the above mentioned physical quantities for the transition. Section [5] consists of similar 
results for the SE-UE crossover. In section [6] we compare our results with semiclassical 
predictions in the weak symmetry breaking regime. Section [7J deals with the phenomenon 
of universal fluctuations, where we compare the variances of conductance and shot-noise 
power for the OE-UE and SE-UE crossovers. We conclude in section [8] with some general 
discussion. Derivation of the skew-orthogonal polynomials and the associated two-point 
kernels are given in the appendices. 

2. The Brownian Motion Model 

In this section we first consider transitions to all three invariant ensembles OE, UE, SE. 
We use the parameter (3 with values 1, 2, 4 respectively for the three cases. In the later 
part of this section as well as in the rest of the paper we specialize to the case of /3 = 2. 

The scattering matrix S for the conductance problem has the standard 
decomposition in terms of reflection matrices r, r' and transmission matrices t, t' as 



Here Ni and N2 refer to the number of modes in the left and right leads connected 



[HE1E2] 




(1) 
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to the cavity. So the matrices r NlxNl and r' N xN correspond to the reflection from 
left to left and right to right respectively. Similarly t N2xNl and t' N xN represent the 
transmission from left to right and right to left respectively. The scattering matrix 
S is thus N s = Ni + N 2 dimensional. Also we write N = min(AT 1 , N 2 ). As a 
consequence of unitarity of S, the Hermitian matrices tH,t'H',l — r^r and 1 — r'V 
have N common eigenvalues T 1; ...,T N with values between and 1. We denote the set 
of these transmission eigenvalues as {T}. 

For considering the transition problem it is assumed that the evolution of S takes 
place as 

S(r + 5r) = S(r)e iv ^ M(r) , (2) 

where M(t), independent for each r, belongs to one of the invariant classes of Gaussian 
ensembles of random matrices [36, 38J. It can be shown [20] that the jpd of transmission 
eigenvalues V({T}; r) obeys the following Fokker-Planck equation: 

«Qpi>=-£P( {r };r). (3) 
In the above equation the Fokker-Planck operator £ is given by 

.7 = 1 J J 



\^ N m)\\{Ti Nl - N ' A+l)/ \l-T k f^\ (5) 



with 

W = -In 

k=l 

Here An({T}) = Y\j <k {Tj — T k ) is the Vandermonde determinant. The equilibrium jpd 
is obtained from CV^ = and is given by [3l [32] 

N p 2 

mm) = Cf\A N ({T})\^]jTf N ^ N2l+l ^. (6) 

3=1 

To make contact with Jacobi weight function we make the transformation 

Xi = 2T 3 - 1, (7) 

and then < Tj < 1 implies —1 < Xj < 1. Thus we get the jpd Pjq ({a?}) for the x 
variables, 

jv 

Pg({*}) = ctf\A N {{x})\e]l{i + x j )^-i\ (8) 

0=1 

where we have introduced 

26 + 1 = |JVi- jV 2 |. (9) 

The Jacobi weight function involved in the above jpd is now easily recognized as 
w n,v{ x ) — (1 — + x) u , with \i = and v = [3{b + 1 — 1//5). In appendix A, 

we give the relevant results for the Jacobi polynomials. The similarity transformation 
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£ = (Peel ) X I 2 P in the Fokker-Planck equation leads to the following Schrodinger 
equation for r) in imaginary time it: 

da{x};r) 



Or 

where 



-HZ({x};t), (10) 



U = (P^y^C(P^ (11) 

is Calogero-Sutherland type of Hamiltonian. Here C is given by (J5J), (JZI) in the Xj 
variables. 

For transitions to UE we have (3 = 2, and this leads to the Jacobi weight function 

w , 2b+1 (x) = (1 + x) 2b+1 (12) 

in (|S]). Moreover in this case the Hamiltonian H splits into N single- fermion 
Hamiltonians minus the ground state energy of the N fermions described by H x , viz. 
U = J2j=i %y - £o- Here 

^-*&-^-H)\tZ) + H)]- (13) 

The eigenfunctions of T-L x are the weighted Jacobi Polynomials WQ^ + i/2(x)P^' 2b+l (x) with 
eigenvalues 

e n = n(n + 2b + 2). (14) 

Also £ = J2n=o £ «' so that the ground state energy of % is zero. 

The formal solution of (E} in terms of the Xj variables is given by 

P({x};r) = e- CT P({x};0). (15) 

Then using the above results, the jpd of eigenvalues for OE-UE and SE-UE transitions 
can be written as [37] 

N 

P({x};T) = C N A N ({x})Pf[F$]l[w(x i ), (16) 

where is the normalization, Pf represents Pfaffian and w(x) is the initial weight 
function. In sections H] and we choose w(x), such that the initial jpd corresponds 
to (jHJ) with P = 1 and 4 respectively. is an antisymmetric function which can be 
obtained from its r = counterpart using the one-body operators introduced in |37J; 
see appendix B. The indices j, k take the values from 1 to N or iV + 1 depending on 
whether N is even or odd respectively. When N is even is simply sgn(xj — Xk)/2 
and — S'(xj — Xk) respectively for the orthogonal and symplectic cases. For odd iV one 
has to introduce additional term Fj T ^ +1 for OE-UE crossover as discussed in appendix 
B. In the SE Kramers degeneracy requires the levels to be doubly degenerate; thus we 
choose both Ni and N 2 to be even so that N is always even. 
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3. Correlation Functions and Conductance Properties 

For the calculation of n-level correlation function i? ra (xi, x n ; r), given by 

N\ f 1 f 1 

R n (xi,...,x n ;T) = 7 — -r ■ I P({x};r)dx n+1 ...dx N , (17) 

we need the two-point kernels (x, y), A^)(x, y) and B^\x, y), which in turn depend 
on the weighted skew-orthogonal polynomials (f)f\x) and the integrated (dual) functions 
ip\ T \x) [37J. These can again be obtained from their initial counterparts (r = results) 
[391 HQ]. See appendix C for the skew-orthogonal polynomials and appendix D for the 
kernels. The polynomials are worked out explicitly in appendices E and F respectively 
for OE-UE and SE-UE transitions, and the corresponding kernels are given in appendices 
G and H. 

The skew orthogonality relation and application of Dyson's theorems [33j E] l ea d to 
the quaternion determinantal expression for R n involving the above mentioned two-point 
kernels; see (ID.5j) . The n-level correlation function lZ n for the transmission eigenvalues 



is related to R n as 

K n {T u .., T n ; r) = 2 n R n (2T\ - 1, 2T n - 1; r) . (18) 

We now consider the application of these results to the quantum conductance 
problem and find expressions for the quantities of interest, viz. the averages and 
variances of both conductance and shot-noise power. The calculation of these quantities 
involve only the level density Ri(x; r) and the two-level correlation function R2(x, y; r). 
These are given in terms of the above kernels as 

R 1 (x;r)=SP(x,x) (19) 

and 

R 2 (x, y; t)=R, (x; t)R 1 (y; r) - S$ (x, y) S% ] (y,x) + A% ] (x, y)B% ) (x, y) . (20) 
We also need the following moments of R\ and R 2 , 



(x) = J xRi(x;r)dx } (21) 
(x 2 ) = [ x 2 R 1 (x]T)dx, (22) 



-i 

(xy) = I / xyR 2 (x,y,T)dxdy, (23) 



-l j -i 

along with the normalizations 



J Ri(x; r)dx = N, 

j j R 2 {x,y;T)dxdy = N{N-l). 



(24) 



(25) 
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These moments along with similar averages (x 4 ) and (x 2 y 2 ) (needed for the variance of 
shot-noise power) can be calculated using i?i(x;r), R 2 (x,y;T) given by f fl9|) , f )20|) and 
the kernels given in appendices G and H. This requires some tedious algebra involving 
repeated use of the orthogonality and the three-term recurrence relations for the Jacobi 
polynomials as given in appendix A. We use these quantities in sections H] and [5] for 
OE-UE and SE-UE crossovers. 

Landauer-Biittiker formalism [H [3J, [HI [151 HI E] enables us to determine a variety 
of transport properties of the mesoscopic cavities using the statistics of transmission 
eigenvalues. Of particular interest are the observables which are linear statistics on 
the transmission eigenvalues, such as conductance and shot noise. These are quantities 
which do not contain products of different eigenvalues. We start with the dimensionless 
conductance which, at zero temperature, is given by the Landauer formula [U HH [15] 



Next we consider the shot noise. These are the time-dependent fluctuations caused due 
to the quantum nature of electron charge, and persist down to zero temperature. The 
expression for the shot-noise power (dimensionless) is known due to Biittiker [31 [161 E] 

as 



N 




(26) 



Thus it follows that the average conductance is given by 




(27) 



N 




(28) 



From this we obtain the average shot-noise power 




(29) 



The variance of conductance is given by 




(xf) (xi) 2 (xxx 2 ) 



(30) 



4 4 4 



Finally we have the variance of shot-noise power, 




(xf) (xl) 2 , (x{xl) 



(31) 



16 16 16 



8 



The correlation functions for OE, UE and SE have been given explicitly in terms 
of Jacobi polynomials in [391 HO]- Using these one can work out the above mentioned 
averages which have been obtained earlier by other methods [U [3j [181 EE] . We find 



N t N 2 



(9) = — -^7i (32) 



{P) (N s -2 + l)(N s -l + l)(N s -l + j j y [66) 

^ /3(iV s - 2 + |)(iV s - 1 + lf(N s - 1 + |) • W 

As pointed out in [19J we find that the exact expression for var(p) does not have a 
compact form for — 1,4 for arbitrary Ni,N 2 . For (3 = 2 however the expression is 
somewhat simpler, given by 

N?N%(N! - l) 2 (N 2 - l) 2 



var(p) = 
+ 



(N s - 3)(N S - 2) 2 (iV s - l) 2 (iV s ) 2 (iV s + 1) 

NfNjjm + 1) 2 (N 2 + l) 2 
(N s - 1)(N B )*(N B + l)2(iV s + 2)2(iV s + 3) 



+ (iV s -2) 2 (iV s -l)(7V s )2(iV s + l)(iV s + 2) 2 - 1 ) 

Also when N\ = N 2 = N, the exact expression for /3 = 1 is given by 
_ N(N + l)(8iV 5 + 6(W 4 + U2N 3 + 91N 2 - A9N - 36) 
Vai ^ ~ 4(2iV - 1)(2N + l) 2 (2iV + 3) 2 (2iV + 5)(2iV + 7) ' ^ ' 

and for /3 = 4 by 

_ N(2N - l)(128iV 5 - 480iV 4 + 568A 3 - 182iV 2 - 49iV + 18) 
Vai ^ ~ 4(4iV - 7)(4iV - 5)(4iV - 3) 2 (4iV - l) 2 (4iV + 1) ' (3?) 

For large Ni, N 2 we can write the compact /3-dependent expression for var(p) correct to 
Oil) as 

2 N 2 N 2 

var(p) = -^[2iV 2 iV 2 + (N, - N 2 )% (38) 

which agrees with the result of [T9] . 

We remark that these relations will also be obtained in sections H] and [5] as r — > 
and t — > oo limits of the transition results. Note moreover that the initial jpd for 
SE-UE transition contains explicitly doubly degenerate levels and therefore the above 
quantities should be suitably modified for (3 = 4 when used in section |5j 



4. OE-UE Crossover 



We consider here the OE-UE transition and calculate the above physical quantities as 
functions of Brownian motion parameter r. The initial jpd of eigenvalues (in terms 
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Figure 1. The density of transmission eigenvalues p(T) = 1Z\(T)/N for OE-UE 
crossover for several values of r for (a)iVi = 1,N% — 1, (b)iVi = 1,N2 = 4, 
(c)Ni = 2,N 2 = 2, (d)JVi = 3, N 2 = 4, (e)iVi = 4, 7V 2 = 20, (f)iVi = 20, 7V 2 = 20. 



of Xj) for this crossover (corresponding to OE) can be found by setting (3 = 1 in (JH]). 
Thus the weight function appearing in (JIB]) in this case is ty(a;) = ty ,fe(x) and the 
crossover in the weight function is from wo t b(x) for OE to wo,2b+i(x) for UE. The jpd of 
eigenvalues for the transition can thus be read from ( Tl6i) . Moreover all the two-point 
kernels can be worked out exactly giving thereby the correlation functions of all orders. 
See appendices D and G for explicit expressions. Using these one can work out (x), 
(x 2 ), (xy) and other moments as mentioned earlier. The corresponding expressions for 
the transmission eigenvalues can be obtained by going back to the T variables as in (J7j) 
and PI). 

The level density 1Z\(T] r) (normalized to N) for transmission eigenvalues is found 
to be 

N-l 

U^T; r) = T 2b+1 + 2b + 2)P°' 2b+1 (2T - 1)P°' 26+1 (2T - 1) 

At=0 
N-l oo 

- T 2b+l + 2b + 2)(-l) N+v e-^- £ » )T P^ 2b+ \2T - l)P°' 2b+1 {2T - 1), (39) 

^=0 u=N 

valid for both even and odd N. Here, on the RHS, the r-independent term, i.e. the 
first term, corresponds to the UE. The r-dependent term is the correction term which 
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contributes in the non-equilibrium regime (0 < r < oo). This term is what leads to the 
weak-localization correction in the average of quantities which form linear statistics on 
transmission eigenvalues. For r — > oo, this term disappears completely and we retrieve 
the well known expression of level density for the UE [301 [31] . 

In the weak-localization study of relevant physical observables in mesoscopic 
systems the quantities are expressed in terms of diffuson and Cooperon contributions 
[22] 123] |2l] . These correspond respectively to the time- reversal symmetry independent 
and dependent parts. The r-independent and dependent terms in ( 139|) can be viewed 
as the analogues of the diffuson and Cooperon contributions. 

For Ni,N 2 » 1 and with fixed ratio N1/N2 the density is given for all r by [3] 



(Ni + N 2 ) T - T c 
K 1 (T) = { 2nT Vl-T' ( 40 ) 



0, otherwise, 



where 

2 



T c = 1^—^ ■ (41) 



N x + N 2/ 

The densities for various values of N\ and N 2 are shown in figure 1 as functions of 
t. In figure 1(e) the densities vanish below some critical T values which are somewhat 
smaller than T c rs 0.44 predicted by (14"!]) ; this is because N± = 4 is not large enough. In 
figure 1(f) the densities for all r are close to that given by (140]) because Ni = N 2 = 20 
is significantly large. 

The average conductance is obtained, by using the value of (x) in (1271) along with 
the substitution N = min(N u N 2 ) and N + 2b+l = max(N u N 2 ), as 

AVVj e _ w , T 

N, N,(N, + 1) K ' 

Recall also that N s = Ni + N 2 . Similarly, substituting the value of (x 2 ) in (129]) gives 
the average shot-noise power as 

{p) = N i N i + n m n ^ - n *? 



(N s - l)(N B )(N a + 1) (N s - 2)(N B )(N B + 1)(N S + 2) 

UV " e-^.(43) 



NMjm - i)(n 2 - 1) o _ 2(Ars - 1)r _ ^(^ + 1)^ + 1) ^_ 2(Ns+1) 



(N s - 2)(N S - 1)(N S )(N S + 1) N S (N S + 1)(N S + 2)(N S + 3) 

The variance of conductance for OE-UE crossover is obtained from (|30]) as 

[9> (N s - l)(N s Y(N a + 1) + (N s - 2)(N S - 1)(N S )(N S + 1) 

NlNj 2NsT ^(^ + 1)^ + 1) 2(Ns+1)T 
(N s y(N s + l) 2 N S (N S + 1)(N S + 2)(N B + 3) 



2N 1 N 2 (N 1 - N 2 



\2 



-N s r 



' (N s -2)(N s y(N s + l)(N s + 2y ' (44) 
The exact crossover result for the variance of shot-noise power can be obtained from 
(13T]) but it is too lengthy to be presented here. We give the asymptotic result in section 

m 
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As expected, these results match with the results (l3"2]) -( l3~4"j) for (3 = 1, 2 in the limits 
t — > and oo respectively. For Ni = N 2 , (f4*2l) and (j4*4j) give back the corresponding 
results of Frahm and Pichard |21j. 

It is clear from the results fj42l - (|44T) that N s t serves as the natural transition 
parameter for the problem. Thus for large Ni,N 2 non-trivial crossover is obtained 
when r ~ 1/N S . The fact that the transition parameter scales as r ~ l/iV s is not 
surprising. It is known, for example, in circular and Gaussian ensembles [361 H2] that 
the global properties such as the level density and variances of low-order traces scale in 
a similar way. This is in contrast to the local fluctuation properties, for example the 
spacing distribution and the number variance, where the transition occurs for r ~ l/N^ 
[34"] 135] |3"6"_[3"?1 |4~2] . We also mention that the r-dependent terms give rise to 0(1/ N s ) 
correctionsjl (the weak-localization corrections) in the respective leading terms of ( 1421) 
and (jl3]). 

The symmetric case of Ni = N 2 » 1 has been studied by many authors. The 
average conductance in this case is given by 

(9) = y - \e~ 2N \ (45) 

which yields the universal value —1/4 for the weak- localization correction in the limit 
t — > [31 [ini [21]. Also in the expression of average shot-noise power given by ( |4"3]) . 
the first r-dependent term vanishes identically, whereas the last two terms cancel each 
other. As a result no weak-localization correction is obtained in the average shot-noise 
power. This is in agreement with the earlier predictions of absence of weak-localization 
correction in average shot-noise power in the symmetric case [31 [TT] . 

The variances of conductance and shot-noise power in the asymptotic limit lead to 
universal results and is considered in section [71 



5. SE-UE crossover 



In this section we consider the SE-UE crossover. Since we have to take into account 
Kramers degeneracy, the initial jpd cannot be obtained by direct substitution of (3 = 4 in 
OH]). Rather we need to consider Ni — > N\/2 and N 2 —> N 2 /2 along with the introduction 
of S functions to take care of degeneracy. Thus there are total of N/2 eigenvalues at 
t = 0, each of them two-fold degenerate. As soon as r ^ the degeneracy is broken 
and we have N distinct eigenvalues. It turns out that the initial weight function for 
this crossover is wq^+i{x). The jpd of eigenvalues in this case is thus given by ( fl~6l) with 
w(x) = wo } b+i{x). Also the expressions for the two-point kernels can again be obtained 
exactly and then R n can be written down (see appendices D and H). 

We get the following expression for the level density 1Z\ (T; r) for transmission 

\ The correction term itself is of 0(1). By C(l/iV s ) correction we mean 1/N S with respect to the 
leading term. 




Figure 2. The density of transmission eigenvalues p(T) = 1Zi(T)/N for SE-UE 
crossover at various r for (a)Ni = 2,N 2 = 2, (b)Ni = 2,N 2 = 4, (c)iVi = 4,N 2 = 4, 
(d)JVi = 4, iV 2 = 8, (e)iVi = 4, iV 2 = 20 and (f)iVi = 20, N 2 = 20. 



eigenvalues in this case: 

2V-1 

T^i (T; r) = T 2b+1 + 2b + 2)P^ 26+1 (2T - 1)P°' 26+1 (2T - 1) 

ju=0 

AT-l oo 

- T 2b+1 ^ ^(2z/ + 26 + 2)(-l) Ar +^ e -( e "-^)^ 26 + 1 (2r - l)P !y °' 2fe+1 (2T - 1). (46) 

^=0 v=N 

Again as in the OE-UE case the first term in the above equation corresponds to the 
UE and the second term is responsible for the weak- ant ilocalizat ion corrections. The 
similarity between (1391) and ( f4"6"l) is noteworthy; the two expressions are identical looking 
except for the interchange of /x and v in the factors (2/x + 26 + 2) and (— l) v inside the 
double summation. For large Ni,N2 and fixed N1/N2, we obtain again ( I40p . 

Figure 2 shows the densities for different Ni,N 2 values as a function of r. Again, 
as in figure 1(e), the densities in figure 2(e) vanish below certain critical values. Also, as 
can be seen from figure 2(f), for large N±, N2 the density is again dominated by the first 
term corresponding to UE. However, in this case the oscillations do not die as rapidly 
as in the OE-UE crossover. 

We now give the results for the average conductance, average shot-noise power and 
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variance of conductance for this crossover. The average conductance turns out to be 

(g) = ^+ NlN2 e-»" (47) 

The average shot-noise power is obtained as 

, v = NfNj N^jN, - N 2 f _ NsT 

W (N s -1) (N B ) (N 8 + l) (N s -2) (N s - 1) (AQ (N s + 2) 6 

- 1)(N 2 - 1) 2(jVs _ 1)T _ Ayv^ + iXiv^ + i) 2(JVs+1)r 

(N s - 3)(N S - 2)(N S - 1)(N S ) (N s - 1)(N S )(N S + 1)(N S + 2) 



(48) 



The variance of conductance comes out as 

, . = N*N* ^(^-1)^-1) (Na _ 1)T 

{y) (N s - l)(iV s )2(iV s + 1) + (N s - 3)(JV B - 2)(JV B - 1)(JV 8 ) 

-2N s r , ^(^ + 1)^+1) 2(Ns + l)T 

(N s - l)2(iV s )2 + ( Nb - l)(N B )(N a + 1)(N B + 2) 



2N 1 N 2 {N 1 -N^ 2 



V -N S T 



(N s -2)(N s -l)(N s )*(N s + 2y • (49) 
Again the exact expression for variance of shot-noise power is too complicated. We give 
the asymptotic result in section [71 

Similar to OE-UE crossover the r-dependent terms lead to 0(1/N S ) corrections in 
the respective leading terms in (1471) and (PI8]) for large Ni,N 2 . The sign of correction, 
however, is opposite to that in OE-UE crossover, hence the term weak antilocalization 
is used. 

In t = limit (]4Tj) -(l49l) agree with the known results (|32l) - (l34l . provided two-fold 
degeneracy for each eigenvalue is properly considered. For that, along with /3 = 4, we 
have to set N x -»■ N 1 /2, and N 2 -> N 2 /2 in (|52" ]l - (fM|l . and then multiply the results for 
(fl 1 )) (p) by 2 and the results for var(p), var(p) by 4 respectively. In the r — > oo limit the 
degeneracy no longer holds so the above results match with (132|) - (!34|) for f3 = 2 without 
any modification. 

We now point out an interesting observation for the OE-UE and SE-UE crossover 
results for (g), (p) and var(g). One can obtain the results for SE-UE from OE-UE results 
and vice versa by simultaneously changing the sign of N±, N 2 and r, and introducing an 
overall — ve sign for (g), (p). For var(g) one does not need to introduce the overall — ve 
sign. 

6. Comparison with semiclassical results 

Weak time-reversal symmetry breaking condition is achieved when r << 1/N B , r = 1/N S 
being the natural crossover scale for the problem. In this case, with Ni,N 2 » 1, we 
get from (T42l) . ()43l) and (T47j) . (pfBl) the following results for the average conductance and 
average shot-noise power for the two transitions: 

(g ^~ T A?(l + iV,T) - (50) 
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N?N$ NiN 2 (Ni — Ni 



\2 



wUr) ■ (51) 

In the above two equations the upper and lower signs correspond to the OE-UE and 
SE-UE transitions respectively. These results should be compared with the semiclassical 
predictions (for the OE-UE transition) for the average conductance and average shot- 
noise power when a weak perpendicular magnetic field is applied to the system [23 ES] , 
viz. 

to> ™_ +c ay (52) 



N, A?(l + C) \K 
N}NI N 1 N 2 (Nj - N 2 ) 2 /1\ 

W = ^ + NHl + C) +0 {nJ' (53) 

where £ is proportional to the square of magnetic flux $ through the system. (We 
have dropped here the extra factor of 2 appearing because of spin degeneracy on the 
RHS in [26J.) We know from the analysis done in |21J that r is proportional to $ 2 
for the quantities considered here. Thus our results are consistent with semiclassical 
predictions. We believe that the semiclassical analysis for the SE-UE transition will 
also agree with the corresponding results in (150|) . fl5Tj) . 

7. Universal fluctuations for conductance and shot-noise power 

We have, from (jUJ) and flUJ) , 

N 2 N 2 

var(,) = ^(l + e- 2 ^), (54) 

s 

for N±,N 2 >> 1. This is valid for both the transitions. Equation ( 154)) leads to the 
phenomenon of universal conductance fluctuations [TJ [21 [TO] in the r = and oo limits 
yielding the universal values 1/8 and 1/16 respectively. Identical result for var(g) for 
(3 = 1 and /3 = 4 arise in the large Ni,N 2 limit because we have explicitly considered 
the Kramers degeneracy in SE. It indicates that the variance of conductance in the 
universal regime does not depend on the spin-rotation symmetry and is sensitive only 
to the time-reversal symmetry. Figure 3 shows comparison of the asymptotic result (1541) 
with the exact OE-UE and SE-UE crossover results. 

We now consider the variance of shot-noise power. As mentioned earlier, the exact 
expression for var(p) for the crossovers is quite complicated and lengthy. However for 
Ni, N2 » 1 it has a simple form. We find 

var(p) = + e — ) + + e — ), (55) 

valid as in ( I54p for both the crossovers. In the r = and 00 limits it leads to the 
universal values 1/64 and 1/128. The exact and asymptotic results for the variance of 
shot-noise power are compared in Figure 4. 

An interesting quantity which involves (g) , (p) and var(g) collectively was suggested 
in [T8]. The ratio NiN 2 veir(g) / (g) (p) assumes the value 2/(3 for all Aq, N 2 . It is natural 
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0.12 



> 




0.08 



0.06 



Figure 3. The variance of conductance as a function of t for various N(= N\ = N2) 
for the two transitions. The dotted horizontal lines represent the universal values 1/8 
and 1/16. 



0.016 



Exact (SE-UE) 



— Asymptotic 




0.04 



Figure 4. The variance of shot-noise power as a function of r for various N(= Ni — 
N 2 ) for the two transitions. The dotted horizontal lines represent the universal values 
1/64 and 1/128. 



to ask the behavior of this quantity for more general cases, for instance for the crossover 
ensembles. It is clear from the expressions of (g) , (p) and var(g) that in the crossover 
regime this quantity has, in general, a complicated dependence on r. However for 
Ni, N 2 » 1, one obtains a rather simple form, similar to above expression for var(g), 

iVxiV 2 ^4 = l + e- 2 ^ (56) 
(9) (P) 
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Note that taking into account the degeneracy for SE leads to identical expressions for 
the two crossovers. Figure 5 exhibits this quantity for several Ni,N 2 values. The OE- 
UE and SE-UE curves have been calculated using the exact results and the dotted curve 
represents the asymptotic results calculated using equation (|56|) . As can be seen this 
asymptotic result works nicely for Ni and N 2 as low as 4. 




— OE-UE 

— SE-UE 
Asymptotic 



W,= 20, 7V 2 = 30 



0.1 0.2 0.3 0.4 0.5 

i 

Figure 5. NiN2va,r(g)/ (g) (p) vs. r for various Ni,N2- The asymptotic curves have 
been drawn using (|56|) . 

8. Conclusion 

In conclusion, we have shown that the transmission eigenvalues in the OE-UE and 
SE-UE crossover regimes are described by the Jacobi crossover ensembles of random 
matrices. We have obtained exact results for the jpd of eigenvalues and the n-level 
correlation functions for these crossovers. By applying these results we have given 
expressions for averages and variances of conductance and shot-noise power as functions 
of Brownian motion parameter r. The relation of this parameter to the magnetic flux 
through the system has been discussed in detail by Frahm and Pichard [21] . 

The exact random-matrix expressions for the averages and variances of shot-noise 
power for arbitrary N\, N 2 have become available only recently [HI [19]. The r-dependent 
expressions derived here not only generalize the results for the averages and variances 
of conductance and shot-noise power to the crossover ensembles but also serve as 
alternative proofs for the — 1,4 and (3 = 2 results in the r — > and r — > oo 
limits respectively. 

Finally we have analyzed the results for large number of modes in the two leads. 
We have obtained the weak-localization corrections in the OE-UE case and the weak- 
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antilocalization corrections in the SE-UE case for the average conductance and average 
shot-noise power. These are in agreement with the corresponding semiclassical results 
in the weak time-reversal regime of OE-UE crossover. 
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Appendix A. Jacobi polynomials 

Jacobi polynomials Pj ,v (x), are defined with respect to the weight function w^ix) = 
(1 — xY l (\ + x) v in the interval — 1 < x < 1 as [13] 

1 w^ v {x)P^ v {x)P£ v {x)dx = hf v 8 jk . (A.l) 
i 

Here hj' v is the normalization and is given by 

k w_ ^ +u+l r(j + / i + i)r(j + ^ + i) 

In the appendices that follow we use the notation 

V n (x) ee P"J b+ \x) (A.3) 
for compactness. V n {x) should not be confused with jpd in (El). Similarly we use 

o26+l 

hj = h 3 = u + b + 1 y (A-4) 

In the main text, however, we stick to the original notation. 

The three-term recurrence relation [13] satisfied by Jacobi polynomials simplifies 
to the following for V n (x): 

xVn[x) ~(2n + 2b + 2)(2n + 2b + 3) n+l[X} 

(26 + 1)2 VJ X) 



(2n + 2b + l)(2n + 2b + 3) 
(2n)(n + 2b+l) 



' (2n + 2b+l)(2n + 2b + 2) Vn - l{x) - (A ' 5) 
This, along with the orthogonality relation flA.lj) . is used in calculating the moments 

d2TD-(ESD. 

Appendix B. One-body operator O x and the antisymmetric function Fjk 

The one-body operator introduced in [37], which is used in finding the expressions 
appropriate to the crossover from r = results, is given by 

O a = V^+^W . w & (B.l) 

W ( X > \/wo,2b+l(x) 
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Here w(x) is same as that in (TTSj) and, as explained earlier, equals Wq^(x) and w 0: b+i(x) 
for OE and SE respectively. Also Wo t2 b+i(x) is the weight function corresponding to 
the UE. Since % x has the eigenfunctions w 0j b+i/2(x)V n (x) with eigenvalues e n , the 
eigenf unctions of O x are (wo^b+iix) / w(x))V n (x) with eigenvalues e~ e " T . Similarly the 
eigenfunctions of the operator (Oj,) -1 are w{x)V n {x) with eigenvalues e £nT . 

The antisymmetric function in f|T6|) can be obtained from its r = counterpart 
by using the operator O. For SE-UE crossover and also for the even iV case of OE-UE 
crossover it equals (xj, Xk) with j, k = 1, 2, N. In terms of the functions dual to 
skew-orthogonal polynomials (see appendix C) it is given by 

oo 

gV{x, y ) = Y&$(*)4>&.i(v) - i>% + i(x)41(y))- (B.2) 

Ai=0 

Q( T \xj,Xk) is obtained by application of Xj and Xk on its r = counterpart, viz., 

g {T \xj,x k ) = Xj O Xh Q m {xj,x k ). (B.3) 
For the iV odd case of OE-UE crossover we also need to introduce 



u y ' (x) = — ' , 



J2(-ire-^V,(x), (B.4) 



then the expression for jpd of eigenvalues (TIB]) holds with the addition F^ +1 = 

Note that for the odd N case of OE co^(x) = 1. To obtain u^(x) from this we 
use 1 = /_j 8{z — x)dz and the expansion 

X 

5(z - x) = wo,b(^)wo,b+i(a;) ^ j-^( z )^i x )- ( B - 5 ) 

At=0 ^ 

Integration over z then leads to ( 1B.4I) with r = 0. Operation of O x on this finally leads 
to the expression valid for arbitrary r (see (1C.6I) ). 



Appendix C. Skew-orthogonal polynomials 

The weighted skew-orthogonal polynomials 0^ (x) and the dual functions ip^ (x) satisfy 
the following skew orthogonality relation: 

j\<f\x)^\x)dx = Z jk , (C.l) 

where Zj k = —Z k j equals 1 for k = j + 1 with j even, — 1 for k = j — 1 with j odd, and 
for \j — k\ 1. Here 

4\x) = j l g^\x,y)^\y)dy. (C.2) 

For the odd A case of OE-UE crossover in addition to the above skew orthogonality 
relation we also have the extra condition, similar to the r = case.[39l HQ] 

j\^(x)<p { ;\x)dx = 5 hN ^ (C.3) 
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because of the presence of the unpaired <f%l.i(x). 

These r-dependent functions are obtained from their r = counterparts [391 HQ] by 
application of the operators O x and (Oj,) -1 as 

# } («) = (oirvf (*), (c.4) 

^ T) (a?) = 0.^f(aj), (C.5) 
wW(i) = 0,w®(i). (C.6) 



< (x, i/) = £ (*)V$ (v) - < } (y)l 



/ N — c \ 



(D.l) 



Appendix D. Two point kernels and the n-level correlation function 

Let 

0, iV even, 

1, N odd. 

The two point kernels are then defined in the following manner: 

jt=0 

+ c&(»)« |r) (y) 1 (D.2) 
a%\x,v)= J2 [^(^^(^-^(^^(y)], 

iu=0 



(D.3) 



+ c - ^Uy)co^(x)}. (D.4) 

Note that for SE-UE crossover N is even only. These kernels can be obtained by applying 
the operators {O x )~ l and O x on the kernels for r = using flC.4j) . (IC.5j) and (IC6|) . Now 
the n-level correlation function can be written in terms of a quaternion determinant as 



i?M(x l5 ...,x n ) = Qdet( j T r^ k { y T ,r j} % i • 

j,k=l....,n 



Appendix E. Skew-orthogonal polynomials for OE-UE crossover 

In this case we have G^°\x,y) = sgn(x — y)/2. The weighted skew-orthogonal 
polynomials <pj(x) and the dual functions ipj(x) for the initial weight function 
wo t b(x) (corresponding to OE) can be compactly written down in terms of the Jacobi 
polynomials Pj' 2b+l (x) using the results given in [39], HO]. The eigenf unctions of O x 
however involve Vj(x). We therefore need to expand (f)f\x) and ipj(x) in terms of 
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Vj(x) prior to the applications of operators to calculate the corresponding expressions 
for the transition. Such expansions are given in [13]: 



1,26+1, 



X 



(l-x)P l n 2b+ \x 



n + 2b + 2) 

2(n + l) 



(2n + 26 + 3) 
Thus, using (1C.4I) and ( 1C.5I) we find, for even iV 



^(/i + 6 + l)P M (x), 

[P„(z) -7Vu(z)] 



2m 



ju=0 ^ 



(r) _ w ,n-i(a;) 

%>mW - 2 5 +1/2 



£ (-l)M e -^ 



/i=2m+l 



^2m+l 



.1' 



_2^+i/2 



-P 2 m(x) + 



2m 



L 2m+1 



-V 2m +i(x) 



2 6+l/2 



3 — E2mT 



V 2m (x) - e 



-£2m + lT 



V 2m +i(x) 



for m — 0, 1, 2, ... . When iV is odd, we have for m = 0, 1, (N — l)/2 — 1, 



2m+l e 



fj.\ x )i 



2m+l 



At=0 



X) 



The unpaired ^-l^) t urns ou t to be 



o£2m + 11~ 



o£2m + 2T 



-2 b+1 / 2 w , b (x) V 2m+1 (x) + - V 2m+2 {x) , (E.9) 



b 2m+l 



b 2m+2 



WQ,b+l{X) 
2&+1/2 



-£2m+lT 



P 2m+ i(a;)-e- £2 -+ 2T P 2m+2 (x 



AT-1 



^ 1 (x)=2S b W^ 1 -^(x) ) 



along with the if) function 



w 0)b+1 (x) 



2 b+i 



X). 



(E.l) 
(E.2) 

(E.3) 
(E.4) 

(E.5) 
(E.6) 

(E.7) 
(E.8) 



(E.10) 



(E.ll) 



(E.12) 



The t = expressions, from which the above results follow, are contained in the above 
equations. 
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Appendix F. Skew-orthogonal polynomials for SE-UE crossover 

The skew-orthogonal polynomials for SE-UE crossover again satisfy the relation (1C.1I) . 
however unlike OE-UE crossover here we have Q^(x, y) = —8'(x — y). As explained in 
section [5j in this case the r = polynomials correspond to the initial weight function 
wo,b+i(x). Since this weight function does not vanish at x — 1, so extra care has to 
be taken in finding the appropriate skew-orthogonal polynomials which satisfy (IC.lj) . 



To calculate the ip^\x) using j^ l Q ( -°^(x,y)(j)^\y)dy, one has to take into account the 
contribution from boundary term also. Also the following relations turn out to be useful 
in arriving at r = expansions of (j) and i/j involving Vj (x) : 

n (1 _ x ) =--J2^ + b)V,- 1 (x) (n>l), (El) 

K 1 '^) = ( t + 2 2&+ 1 l) [M*)-T>*-i{*)] (»>0), (F.2) 
along with P_"{x) = 0. It is to be noted here that the polynomials P^ 1 ' +1 (x) form 



orthogonal set for \i > 1 only. These relations follow from the results given in [43] . 

Once the r = results are known the results for arbitrary r can be written down 
easily using the operators O x and (Oj,) -1 . The final results for the and ip functions 
for this crossover read 

I n 2m 

= £(-1)^7^), (F.3) 

At=0 

oo —er 

Vfi(x) = -2 b+ V 2 w 0jb (x) 6 — V ^ (F-4) 

H=2m+l V 



,(r) , v _ WQ,b+l(x) 
J 2m+l\ x ) 2 b+1 l 2 



V 2m {x) - e e ^ T V 2m+l {x) , (F.5) 



4 T m+ i(*)= - 2 b+1/2 WoA*) ~ h + - h V 2m+1 (x) . (F.6) 

Appendix G. Two-point kernels for OE-UE crossover 

We substitute the expressions of <fij(x) and ipj(x) in (|D.2j) - flD.4[) and f )B.2[) . By 
manipulation of the summations we obtain the following expressions for the kernels 
which are valid for both even and odd N: 

sP(x, y) = w , b (x)w , b+1 (y) \£ j-V^VM £ S h ^)V u (y) 

(G.l 
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^=1 u=0 



N-l M-l 226+l e (e M +e„)r 



hnK 



oo oo 



b n (x,y) = -w Qib+1 (x)w ,b+i(y)^2 Yl 

H=Nv=fi+l 



oo oo 



Q (r) (x,y) = w 0tb+1 (x)w 0)b+1 (y) 

/i=0 u=fj,+l 
oo /i— 1 

= -w ,6+i(z)wo,&+i(y) 

^=1 i/=0 





2 2b+l 


(-I)' 1 






226+1 


(-l)M- 





V ll (x)VM-V v {x)V ll (y) j 

(G.2) 

V ll {x)V v (y)-V v (x)V lt (y) 

(G.3) 



2 2b+l 



(G.4) 

The first expansion of ^^(x, £/) in (IG.4I) is obtained by using the <fi and ip belonging to 
the even- A" case in (IB. 21) . whereas the second one follows by using <p and ip belonging 
to the odd-A" case. As expected, these two representations for Q^ T \x,y) are equivalent 
and follow from each other by changing the order of summations. 



Appendix H. Two-point kernels for SE-UE crossover 

(-l)"e 



For the SE-UE transition we obtain 

N-l 

S N~\ x ,y) = w 0tb+ i(x)w 0>b (y) ^2 —V^V^y) - ^ 

fi=0 ^ n=0 u=N 



iV - 1 00 ' -> e -( £ „- £M )r 
//„'' / '" : " ^ ^. ~h„ 



(H.l) 



a n( x >v) = w o,b+i( x )wo,b+i(y) Y Yl 



N-l fi-1 



2^+i 



V^x)V u (y)-V u (x)V^y) 



oo oo 



b nH x >v) = W 0,b( X ) W 0,b(y) Y Y2 

H=N v=n+l 



226+l e -(e fl +e„)T - 



V^x)V„(y)-V u (x)V^y) 



(H.2) 



(H.3) 



oo oo 



G {r) (x,y) = -w 0jb (x)w Qtb (y) ^ ^ 

/i=0 i/=fi+l 



226+l g -(e M +£^)r 



VJ x )V v (y)-V v {x)VJy) 



V^x)V v {y)-V v {x)V^y) 



(H.4) 

In flH.4p the first expansion follows from the substitution of 4> and r/> for SE-UE crossover 
in (1B.2j) . We have given the second expansion as an analogy to the corresponding 
expansion in (1Q.4|) by changing the order of summation in the first one. 
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